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Abstract 

We study orbital magnetism of a degenerate electron gas in a number of 
two-dimensional integrable systems, within linear response theory. There are 
three relevant energy scales: typical level spacing A, the energy T, related 
to the inverse time of flight across the system, and the Fermi energy ep. 
Correspondingly, there are three distinct temperature regimes: microscopic 
(T <C A), mesoscopic (A < T < T) and macroscopic (T <C T <C Ep). 
In the first two regimes there are large finite-size effects in the magnetic 
susceptibility x, whereas in the third regime x approaches its macroscopic 
value. In some cases, such as a quasi-one-dimensional strip or a harmonic 
confining potential, it is possible to obtain analytic expressions for \ m t ne 
entire temperature range. 
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1. Introduction 

A degenerate electron gas, in the presence of a weak magnetic field, ex- 
hibits weak orbital magnetism 1 (the Landau diamagnetism). For a two- 
dimensional gas the value of the orbital magnetic susceptibility is given by 
Xl = — e 2 /12nMc 2 , where M is the electron mass. (Double degeneracy with 
respect to spin is implied in this expression, as well as in all subsequent 
formulae.) This result applies to a macroscopic system. 

Whether a sample of a given size L can be considered as macroscopic, 
depends on the temperature T. Namely, T (in energy units) should be 
compared with the characteristic size-dependent energies such as the mean 
level spacing, A ~ 2irh 2 /ML 2 , or the inverse time of flight across the sam- 
ple, T ~ hv F /L ~ k F LA/2-K, where v F and k F are the Fermi velocity and 
wave number. Generally, one should distinguish between three temperature 
regimes: 

(i) For T < A (the "microscopic" regime) discreteness of the energy levels 
comes into play and the sample can be viewed as a giant atom. The 
magnetic response in this case can be very strong and includes such 
exotic possibilities as perfect diamagnetism and Meissner effect 2 . 

(ii) For higher temperatures, A < T<r, the system enters the mesoscopic 
regime (for recent reviews see Refs 3,4). Here the typical value of the 
magnetic susceptibility is of order (kpL) a | \l | and can have either 
sign. The exponent a depends on the sample geometry. For most 
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two-dimensional integrable systems a = 3/2, although other values are 
also possible in some special cases (see below). For completely chaotic 
two-dimensional systems a — 1. 

(iii) For still higher temperatures, when T ^> T (but smaller than the 
Fermi energy ep), the system can be considered as macroscopic and its 
magnetic susceptibility, up to small corrections, is given by the Landau 
value xl- 

Thus, at present there is a good qualitative understanding of the phe- 
nomenon of orbital magnetism in various temperature regimes and for var- 
ious geometries. However, reliable quantitative results are scarce. Most of 
such results refer to the mesoscopic regime 4-6 and are based on a semiclas- 
sical approximation for the density of states. This approximation becomes 
inadequate both at very low temperatures, T < A, and at high tempera- 
tures, T > T. It, therefore, seems useful to consider a few simple systems, 
starting with an exact expression for the susceptibility Xi an d to observe the 
behaviour of x in the entire temperature range. This is done in the present 
paper by using a linear response expression for x- 

In Section 2 we present several equivalent expressions for the orbital mag- 
netic susceptibility within the linear response theory. In Sections 3, 4, 5 we 
consider specific examples of a strip, disc, and square-geometries. Section 6 is 
devoted to an electron gas confined by a two-dimensional parabolic potential. 
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2. Linear Response Theory for the Magnetic Susceptibility 



We consider an electron gas, confined to some domain in the xy-plane and 
subjected to a weak magnetic field B in the z-direction. The grand-canonical 
potential 

n = -I f dEp(E) ln[l + e^-^] , (1) 
p J 

where \x is the chemical potential, (3 = 1/T and p(E) is the single-particle 
density of states. Sometimes it is useful, by integrating by parts twice, to 
rewrite Eq. (1) as: 

n = -fdEQ(E)^ (2 ) 
where f(E) = {1 + exp[/3(E — /i)]} -1 is the Fermi function and the quantity 
Q(E) = - f E dE' f E dE"p(E") (3) 

J— oo J — oo 

has the meaning of a grand-canonical potential, for the same system, at zero 
temperature and with the chemical potential equal to E. 
The density of states can be written as 

p(E) = ——ImTrG(E) = --ImTr (E + irj- H)' 1 , (4) 



7T 7T 



where G(E) is the retarded Green's function, at energy E. The full (single- 
particle) Hamiltonian H is split into the unperturbed part, H a = p 2 /2M, 
and the perturbation 
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which describes the effect of the (static) magnetic field. It is assumed that 
the vector potential A satisfies the condition divA = 0. 

Expanding G(E) in terms of the unperturbed Green's functions, G (E) = 
(E + 17] — H )~ l , one obtains, up to second order, 

G = G + G VG + G VG VG , (6) 

which, on substitution into Eq. (4), leads to the following expression for the 
correction Sp(E) to the density of states: 

S P = ~ (^) 2 ^ImTr{[ Go (A-P)} 2 + M(G A 2 )} . (7) 

The first and the second term describe, respectively, the para- and diamag- 
netic corrections. Plugging Eq. (7) into (1) and integrating by parts gives: 



1 / e x 2 



J dEf{E)ImTr\[G (A-p)] 2 + M(G A 2 )} . (8) 



27T VMCy 

This expression is proportional to B 2 , so that the susceptibility, in the B — > 
limit, is x — —25Vl/B 2 A, where A is the sample area. 

For further calculations it is useful to have an expression for x m terms 
of Xo(E) which is the susceptibility of the system at T=0, ix — E (compare 
with Eq. (2)): 

X = -JdE Xo (E)^ (9) 

Eq. (8) is written in an abstract operator form. For practical calculations 
it is useful to use a particular representation. For example, if one computes 
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the traces using as a basis the eigenstates of H Q (with the appropriate bound- 
ary conditions), one obtains, after some algebra: 

Here the summation is over all states of the unperturbed Hamiltonian H , 
with eigenenergies The energies e[ and e" denote the first (i.e. propor- 
tional to B) and the second (proportional to B 2 ) corrections to the unper- 
turbed energies £j. The first correction can exist only for levels which are 
degenerate in the absence of B. 

Another useful expression for \ is obtained by writing Eq. (8) in the 
coordinate representation. This results in: 



X = 



nAB 2 \McJ 

x l-n 



Im J dEf(E) J d d r J d d r' x 



A(r)~G (r,f';E) 



Mr')--^G (r',f;E 



+ 



+M5(f-f')G (f,f';E)A 2 (f) } , (11) 

where G (f : r';E) is the unperturbed (retarded) Green's function in the co- 
ordinate representation. Since this representation is usually the most ap- 
propriate for making various approximations, Eq. (11) is a good starting 
point in many cases. It was used, for instance, in the study of mesoscopic 
effects in disordered systems 7 (in this case Go includes the random poten- 
tial of impurities). It also helps to prove in the most direct way that, for 
T > T, the susceptibility approaches its macroscopic value xl, indepen- 
dently of the sample geometry. Indeed, the Fermi function f(E) has poles at 
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values E n = fi + i(n/ (3)(2n + l). Therefore the integral over E in Eq. (11) can 
be replaced by a sum containing Go(r,r'; E n ). This function in an infinite 
system decays with distance as exp(— k F \ f— r'\ ir(2n + 1) //?//). It is there- 
fore clear that for a system with size L larger than (3fx/k F , i.e. for T > T, the 
susceptibility ceases to depend on sample size or on its geometry. Therefore, 
for T > T linear response theory is applicable as long as the cyclotron energy 
%uj c is smaller than T. However, for T < T the susceptibility x does depend 
on sample size and its geometry, and the linear response condition requires 
that hu c is smaller than the level spacing A (i.e. the magnetic flux $ through 
the sample is smaller than the flux quantum <l> = 2nhc/e). 

A useful approximate expression for x is obtained upon using the semi- 
classical approximation 8 for the Green's functions in Eq. (11). Let us briefly 
outline the derivation (details are presented elsewhere 9 ). First, one derives a 
semi-classical approximation for the function xo(E). This is done by rewrit- 
ing the Green's functions in terms of the propagators K(f, f', t), approximat- 
ing the propagators by their semi-classical expressions 8 and performing the 
integrals within a saddle-point approximation. Then one substitutes Xo(E) 
into Eq. (9) and integrates over E, making use of the approximation 



where F(E) is some smooth function of E (i.e. does not change appreciably 
within an interval of order T) and 1Z(x) = xj sinhx. The resulting semi- 
classical expression for x is 




(12) 




(13) 
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Here A labels families of primitive periodic orbits for integrable systems or 
isolated orbits for chaotic ones, r is the winding number, t\ is the period of 
a primitive periodic orbit and r t = H/ttT. Factors rf Ar (/i) are related to the 
oscillating part of the (unperturbed) semi-classical density of states 8 ' 10 : 



For integrable system (A\) is an orbit area squared and averaged over the 
family A. For a chaotic system it is simply the squared area of an isolated 
orbit. Averaging over a family amounts to integration over one of the angle 
variables 4 



The semi-classical expression for x, Eq. (13), coincides with the one derived 
in Ref. 4, in the B — > limit. This fact demonstrates that it does not matter 
which of the two approximations, i.e. linear response or semi-classics, is done 
first. 

3. Strip Geometry 

In this section we consider electrons confined to a strip —L x /2 < x < 
L x /2, —Ly/2 < y < Ly/2. Periodic and zero boundary conditions are as- 
sumed along x and y directions respectively, and the limit L x — > oo is taken. 
It is convenient to choose the Landau gauge: A x = —By. A y = A z = 0. Sta- 
tionary states are labelled by two quantum numbers, — oo < k < +oo and 



Posc(/i) =^2d x ,M ■ 



(14) 




(15) 
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n 



1,2,.... The eigenfunctions ip n ,k(x, y) = e u nj k(y), where u n ^ satisfies: 



u n ,k{y) = eu n ,k(y) 



Treating the magnetic field as a perturbation, one obtains 11 that the first 



order correction e' nk = and the second order correction 



"nk 



L y 



eB\- 

2 1 \ r \~) 
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7r 2 n 2 



+ 



2 r2 

2/ 



7T 



15 ^ 

-,4 



(16) 



7r* y rr rr 

Thus, only the first term in Eq. (10) is present, and the zero-temperature 
susceptibility (including spin degeneracy) Xo(E) is given by 

4 r+°° dk 



4 r+co ^ w 



(17) 



where e n k = (h 2 k 2 /2M) + (h 2 n 2 n 2 /2ML 2 y ) are the unperturbed energy levels. 
Next, we integrate over k, apply the Poisson summation formula 1 to the sum 
over n, and insert the resulting expression for Xo(E) into Eq. (9). The integral 
over E is then performed, using the approximation (12), which amounts to 
neglecting small terms of order T / \i. The final expression for x i s: 



X_ 

XL 



1 + w /^g COS ( 2ikFL v - j) n ( ZntT ^ + 



71 



e=i 



£3/2 



+ 0(l/Jk F L v ) , (18) 



where k F = y2Mfj,/h 2 , T = h 2 kF / 'ML y = hv F / 'L y and the function lZ(x) has 
been defined above. In Eq. (18) we wrote down only the leading oscillating 



term, of order ^Jk F L y , and the term 1, responsible for the Landau diamag- 
netism. There exist also small oscillating corrections, of order {k F L y )~ 1 / 2 
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and smaller which are not written down explicitly, even though they are 
calculable. 9 Let us only mention that, in addition to oscillating corrections, 
there is also a small non-oscillating paramagnetic correction, 9\xL\/^kpLy, 
to the Landau value xl- 

Thus, the strip geometry provides a rare example for which it is possible to 
obtain an essentially exact (up to small corrections of order T j \x) expression 
for the linear susceptibility x, including all size-dependent terms. One can 
observe the change of x with T in the entire temperature range, from T = 
and up to T ^> T when x becomes equal to its macroscopic value xl- Eq. (18) 
is quite similar to the corresponding expression for the case of a parabolic 
confinement. 12 This fact demonstrates that the nature of confinement (i.e. 
hard walls or "soft" confinement) is immaterial for the phenomenon of size- 
dependent fluctuations. 

In fact, the leading oscillating term in x can be obtained, within a semi- 
classical approximation, for any confining potential and for arbitrary mag- 
netic field. 9 For small fields, the oscillations are given by an expression similar 
to Eq. (18), up to an overall factor of order unity and an extra phase in the 
argument of the cosine. L y should be understood as some effective width of 
the confining potential. 
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4. Circular Geometry 

The electron gas is confined to a disc of radius R. The unperturbed, 
i.e. zero-field stationary states are given, up to a normalisation factor, by 
exp(im9) J m (X mn r/R), where \ mn is the n-th zero of a Bessel function J m {x). 
The unperturbed energies are e mn = (h /2MR 2 )X 2 nn . A pair of states \m,n) 
and | — m, n) have the same energy. 

The perturbation term in the Hamiltonian, due to the magnetic field, is: 

ifteB_d_ e 2 B 2 2 
and the first and second-order energy corrections are: 

e ' mn = ~Yu\ m ' e ™ n = !^( mn i r2 i mn > • ( 20 ) 

Note that the first-order term in V does not contribute to the correction d' mn 
which is therefore purely diamagnetic. It now follows from Eq. (10) that 

i^S = — 12 (E( mn l r V^)/(w) + T7E m2 ^-| • ( 21 ) 

\XL\ - rt (mn lvl mn ut mn ) 

We analyse first the low-temperature regime, T < A e 2h 2 /MR 2 . Since 
{mn\r 2 \mn} ~ R 2 , the first (diamagnetic) term is of order of the total number 
of electrons, iV ~ (k F R) 2 = 4/x/A. The second (paramagnetic) term exhibits 
sharp peaks every time when the chemical potential /i coincides with an 
energy level e mn . Indeed, for /i = e mn , the function df/de mn is equal to 
— 1/4T, and it rapidly decreases when /z deviates from e mn by a few T. Since 
the quantum number m, for states near /i, is typically of order kpR, the 
height of the paramagnetic peaks is of order (k F R) 2 A/T ~ fi/T. Thus, in 
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the low-temperature regime the susceptibility x (normalised to the Landau 
value \xl\) plotted as a function of //, displays a diamagnetic background, of 
order ///A, with sharp paramagnetic peaks of height fi/T and width T. An 
exact numerical computation of the expression (21) confirms this qualitative 
picture (Fig. 1). 

Next, we consider temperatures in the range AcT^Te Tivf/^R- For 
such temperatures the paramagnetic peaks get smeared and the diamagnetic 
background, of order (kpR) 2 , cancels with the corresponding paramagnetic 
term. The net effect is an oscillating term, of order (/c^i?) 3 / 2 . The calculation 
is based on a Poisson summation of the double sum in Eq. (21) and on a semi- 
classical approximation for the unperturbed energies, or \ mn , which satisfy 13 : 

yj \ 2 mn -m 2 -m arccos(m/ A mn ) = ir(n + -) (22) 

Let us outline the calculation of the paramagnetic term 
X^/\Xl\ = —{Q^ 2 /MR 2 )Y^m 2 (df/de mn ). As usual, it is simpler to con- 
sider first the zero-temperature case and then to use Eq. (9). At T = 0, 
df/de— —5(fi — e) and 

X? ) (^) = 12| X l|E^L-7 2 ), (23) 

m,n 

where 7 2 = 4/i/A. After applying the Poisson summation formula, m and n 
go over into continuous variables: m — > x, (n + |) — > y. The integral over 
y is immediate, due to the 5-function. The integral over x is done in the 
saddle-point approximation, using the large parameter 7 and the continuous 
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version of Eq. (22). The resulting expression for Xo\E) is: 

X^\E) = \xl\1\i 2 + 1 V2 E <KK X , K v ) + Ofr)) , (24) 

l 4 K X ,Ky j 

where the sum runs over 1 < K y < oo and < 2K X < K y , and 



^Ky 



COS 



2 7 ^sin^j-^ + ^ 



(25) 



A similar treatment of the diamagnetic term results in a term — 37 2 /4 plus 
oscillating corrections of order yfiy. Thus, the large non-oscillating terms 
cancel and the net result for Xo(E) is given by the second term in Eq. (24). 
Finally, using Eq. (9), (12), we find: 



x = \xl\i 3/2 E <f>(K x ,K y )1l( 

K X ,Ky \ 



n KyT . / K x \\ 
-^sm^-jj. (26) 

Thus, in the temperature range A < T < T, the susceptibility x 1S , typically, 
of order (kpR) 3 ^ 2 and can have either sign. It oscillates, as a function of /x, 

1 /2 

with a period of order (/iA) ' ~ T. In Fig. 2 we present x as a function of 
7, as obtained from Eq. (26) (solid line). For comparison, dots represent the 
result of a numerical computation, based on the exact Eq. (10), with e- and 
e" given by Eq. (20). These oscillations reflect the structure of the density of 
states smoothed over energy intervals AE < T, in the same way as the sharp 
peaks of the low-temperature regime reflected the exact (discrete) spectrum 
of the system. Eq. (26) can be obtained from the corresponding expression 
of Ref. 4, in the B — > limit. The derivation in Ref. 4 was based on a 
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semi-classical approximation for the density of states 8,10 , with a subsequent 
introduction of the magnetic field via the classical action. In contrast, we 
have started with an exact expression for the linear susceptibility and used 
(in the mesoscopic temperature regime) a semi-classical approximation for 
the energy levels of the system. 13 



5. Square Geometry 

Here we discuss electrons within a square of size L. The unperturbed 
energies are e nm = (ir 2 h 2 /2ML 2 )(n 2 + m 2 ). A state \n,m) is degenerate with 
the state \m,n) (there can be, in addition, accidental degeneracies if a pair 
n', m! has the same sum of squares as the pair n, m). 

Let us first consider low temperatures, T<Ae 2n% 2 /ML 2 , and discuss 
the paramagnetic peaks due to the second term in Eq. (10). The first order 
correction e' nm is due to lifting of the double degeneracies by the magnetic 
field. (We do not consider accidental degeneracies, although the treatment is 
readily extended to include that case as well.) The degeneracy is lifted only 
if n and m have different parity and in that case 

, 32 heB n 2 m 2 

turn = -tmn = ^ M C (n 2 - m 2 f ^ ' 

The largest corrections occur when n — m ± 1. In such cases e' nm ~ 
(heB /Mc)kpL and the height of the corresponding peak in susceptibility 
is Xmax — \XL\(k F L) 2 A/T, just as in the case of the disc. Note, however, 
that in the square geometry such large peaks are much more rate than in a 
disc. For a disc large peaks were separated by a distance of order A. For a 
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square such peaks occur, roughly, for each pair (n, n + 1), i.e. are separated 
by a distance of order nA ~ /c^LA ~ T = Hvf/L. The difference between a 
square and a disc is clearly seen, if one compares Fig. 3 to Fig. 1. Except for 
the large paramagnetic peaks in Fig. 3 one can see an oscillating background. 
This background comes form the first term in Eq. (10). It will be shown be- 
low that this term can be either para- or diamagnetic and its typical value 
is of the order (/c^L) 3 / 2 . 

In the mesoscopic temperature regime, A < T < T, the susceptibility x 
is described by the semi-classical expression (13). This case was studied in 
detail in Ref. 6 and particularly in Ref. 4 where expressions for the factors 
d\ tr and (A\) can be found. The resulting expression for x i s: 



where u x and u y are positive coprime integers, which label primitive orbits. 
Only odd u x and u y enter the sum in Eq. (28), since otherwise the area 
enclosed by an orbit is exactly zero. 4 

As an example, in Fig. (4) we present x/\Xl\ as calculated from Eq. (28) 
(dots). We chose the same ratio T/A = 5 as in Ref. 4. The result is practi- 
cally indistinguishable from the numerical one (solid line), based on Eq. (10). 
(Let us note, that in Ref. 4 some disagreement between Eq. (28) and numerics 
was observed). There are some qualitative differences between the mesoscopic 



X 



5v^ 
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(28) 



15 



oscillations in the square geometry (Fig. 4) and those for a circle (Fig. 2): in 
the circle oscillations are modulated on an energy scale which exceeds T by 
an order of magnitude. 

Comparison with an exact numerical computation demonstrates that ex- 
pression (28) is valid down to temperatures T m A, but fails for T < A. 
Nevertheless, it can be used for estimating the aforementioned background, 
due to the first term in Eq. (10). The point is that this term ceases to change 
when temperature is lowered from T ss A down to T = 0. To obtain an es- 
timate of expression (28) at low temperatures, we square it and average out 
the fast oscillations. This gives, for the typical value of x 2 in the background 



The sum over repetitions r is estimated by replacing it with an integral, with 
an upper cutoff provided by the function 71. This gives a logarithmic factor, 
so that 



Let us mention that a similar logarithmic factor appears in chaotic billiards, 
at low temperatures. 5 The behaviour of x i n that case is, of course, quite 
different form the square geometry: the enhancement factor is k F L, instead 
of (kpL) 3 / 2 , and the large paramagnetic peaks disappear. 




(29) 




(30) 
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6. Harmonic Confinement 

Our last example is a degenerate electron gas confined by a harmonic 
potential. 14-17 The problem of orbital magnetism in this case was considered 
previously by a number of authors, and some analytical and numerical re- 
sults have been obtained for various temperatures and fields. Below we shall 
derive an analytical expression for the linear susceptibility Xi as a function 
of temperature. 18 This expression demonstrates the crossover from strong 
magnetic effects towards the weak Landau diamagnetism, under increase of 
temperature. 

We consider a slightly anisotropic harmonic potential, 
U(x,y) = (M/2)(Qix 2 + Q 2 y 2 ), with Qi close to fl 2 , namely 

— Q 2 \ = Af2 < hfll/pL. In the absence of a magnetic field the spectrum 

is given by e nm = Mli{n + \) + %VL 2 {m + |), where n, m — 1, 2, For the 

isotropic energy levels can be labelled by a single integer 

t = 1,2,..., and the £-th level is Mold degenerate. A small anisotropy aQ 
splits the degenerate levels into narrow "multiplets" . The width of the n-th 
multiplet is of order £h,AQ, which is fiAfl/Q for multiplets near the energy 
/i. The above formulated condition, aQ < TiVL 2 / fi, is just the requirement for 
having well defined multiplets. It is clear on physical ground, and is verified 
by the calculation below, that such weak anisotropy can affect the physical 
properties of the system only at temperatures, T < fiAfl/Q. 

When a weak magnetic field is applied, energy levels acquire a second 
order correction 




(31) 
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Since for any finite anisotropy the first order correction is zero, Eq. (10) gives: 

f(e nm ) , (32) 



X 24ft 1 ^ 



|xl| Mi? 2 ft? - ft 2 . 
The effective radius i? is determined from MVL 2 R 2 /2 = //, where ft = (fix + 
Q 2 )/2 ~ fti. Applying to the double sum in Eq. (32) the Poisson summation 
formula, we obtain for the zero-temperature susceptibility Xo(E): 

where A = f^/fti and 7 = /i/ftft = k F R/2. 

Calculating the elementary integrals and making use of Eqs. (9), (12), one 
can obtain a final expression for the susceptibility x(T). This expression is 
rather cumbersome and will not be given here. 9 For temperatures T 3> 7 ft Aft 
and anisotropy Aft <C ft/7 it simplifies to: 

-2L = -l + 27 2 f]cos(27rr7)^ f J (34) 

IXil r =l V 1 / 

where T = Tivf/2R = ftft/2. Except for the leading oscillating term, of order 
7 2 , there are smaller oscillating terms which are not written in Eq. (34). This 
equation does not contain anisotropy, which demonstrates that, for T much 
larger than the multiplet width 7/iAft, the anisotropy does not come into 
play (up to exponentially small corrections). For T ^> T, the oscillations in 
Eq. (34) are negligible and x = Xl- for T < T, one can keep only the first 
term in he sum, which results in an oscillating term of order 7 2 ~ (kpR) 2 - 
For ryhscriptstyleAtt < T <C T, many terms contribute to the sum. The 
resulting expression exhibits paramagnetic peaks, of height {kpR) 2 Y /16T 
and width T, on a diamagnetic background of order —{kpR) 2 - 
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For T < 7^,aQ, Eq. (34) ceases to be applicable. An analysis of the full 
expression 9 for x(T) shows that it matches the Eq. (34) at T ^Tia.VL and 
that only minor changes in x{T) occur when T is lowered below 7ft a£). This 
means that at low temperatures, T < jfiAQ,, the width of the paramagnetic 
peaks becomes ^TiaVL and their height is of order k F RT /HaQ,. 

This result is a manifestation of a general rule. If there is a cluster of 
nearly degenerate levels about some energy e c (the width of the cluster 5 is 
much smaller than the typical level spacing A), then for T > 5 the cluster 
behaves as a single degenerate level: it gives rise to a paramagnetic peak 
of height g(k F L) 2 A/T, where g is the number of non-zero eigenvalues of 
the matrix (i\L z \j) in the subspace of the cluster. For T < 5, the width 
of the cluster becomes relevant and the peak saturates at a value of order 
g(k F L) 2 A/S. This is best seen if e- and e" in Eq. (10) are written explicitly 
using the symmetric gauge. This leads to: 



When 5 < T < A the paramagnetic contribution of the cluster, as given by 
the second term in Eq. (35), is approximately — (3n / AM) f (e c )trL 2 z , where 
the trace is within the subspace of the cluster. Estimating the trace as 




(35) 



gh 2 (k F L) 
peak. 



2 leads to the expression g(k F L) 2 A/T for the Eq. paramagnetic 
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Let us, finally, mention that when the anisotropy Af2 approaches the value 
fi/7 (i.e. neighbouring multiplets start to overlap), the well pronounced para- 
magnetic peaks disappear, although the oscillations are still of order (k F R) 2 . 

7. Conclusions 

We have calculated the linear magnetic susceptibility xCO f° r several two- 
dimensional integrable systems. Generally, there are three distinct energy 
scales: level spacing A, inverse time of flight T and the Fermi energy €f- 

For T < A, the susceptibility is sensitive to the detailed structure of the 
energy spectrum as well as the matrix elements of the angular momentum 
operator. In particular, degeneracies lead to large paramagnetic peaks of 
order g(kpL) 2 A/T, g being the level degeneracy. At such low temperatures 
the sample dimensionality is of no importance and similar peaks exist also 
in three-dimensional systems with degenerate levels. 19 For a pair of nearly 
degenerate levels, when the level separation (5 < A, the height of the corre- 
sponding peak saturates at T m S. The typical value of x between the peaks 
is also non- universal: for systems with rotational symmetry it is of order 
XL{k F L) 2 whereas for a square it is of order \xL\(k F L) 3 ^ 2 ^J\n(k F L) and can 
be either para- or diamagnetic. 

For the mesoscopic temperatures, A < T < T, the orbital magnetic sus- 
ceptibility remains large. In units of \xl\, it is of order (k F L) a and oscillates, 
as a function of the chemical potential /i, with a period a/j, ~ T. For generic 
integrable systems a — | (compare to a = 1 for chaotic systems 3 ). For 
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quasi-one-dimensional systems (a strip) a — \ and for a harmonic confin- 
ing potential a = 2. Thus, the harmonic potential is a very special case 
even among the integrable systems. The point is that, for the isotropic case 
and at zero magnetic field, the two-dimensional harmonic potential exhibits 
"accidental" degeneracies, which are not related to the rotational symmetry 
(alike the "accidental" degeneracies in the hydrogen atom). 

Finally, for T <C T <C //, all large orbital magnetic effects disappear and 
X becomes equal to the Landau value xl- Thus, for the macroscopic limit to 
be achieved, it is not sufficient to have T 3> A, as one might naively expect. 
A much more stringent condition, T ^> T, is needed. We close the paper 
with the following remarks: 

(i) Although our calculations have been done for the grand-canonical en- 
semble, one can immediately infer about the picture for the canonical 
ensemble. For T<A, the susceptibility can be very sensitive to the 
exact number of particles in the system, provided that the last oc- 
cupied level is degenerate. If such a level is partially occupied, the 
system is paramagnetic with x — c \XL\(kpL) 2 A/T, where the factor 
c accounts for the degeneracy and occupancy of the level. For a fully 
occupied level, and in the presence of rotational symmetry, the re- 
sponse is diamagnetic, with x — XL(k F L) 2 . (The discussion can be 
generalised to the case of a cluster of nearly degenerate levels, as was 
done above for the grand-canonical ensemble.) For mesoscopic temper- 
atures, A C T < T, the standard transition from the grand-canonical 
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ensemble to the canonical one applies, i.e. Xam(N) = Xgr(M^O)- I n ~ 



tides is much smaller than the change in the number of particles, 
(aN) ~ ( A/i) (dN/d[i) ~ T/A, corresponding to a significant change 
in x- 

(ii) So far we have discussed only the orbital magnetic susceptibility. The 
Zeeman splitting, due to the electron spin, also contributes. Within 
the linear response its contribution Xs is simply added to the orbital 
part of the susceptibility. This is clearly seen form Eq. (10), if the 
Zeeman splitting ±/ie-B is included into the first order energy correction 
e[. This gives Xs = A~ l fi 2 B {dN/dfi), where N = Eif(^i)- At low 
temperatures, T <C A, Xs exhibits paramagnetic peaks Xs — \xl \ A/T, 
which is just the Curie paramagnetism due to the last occupied level. 
For temperatures T ^> A, Xs is given by the Pauli paramagnetism, 
Xs = 3|xl|, up to small corrections. Thus, in this case, there is no 
appreciable mesoscopic effects due to the electron spin. 

(iii) The magnetic field in this paper was treated as a given (homoge- 
neous) external field, B cxt . For sufficiently large Xi however, the or- 
bital magnetic currents flowing in the sample produce a field -B cmT 
which is comparable to -B ex t- In three-dimensional sample this hap- 
pens when the magnetisation M becomes comparable to _B cxt /47r, or 
|x| ~ I/An. In two dimensions the magnetisation is defined as the 
magnetic moment per unit area (rather than per unit volume), so that 



deed, the thermal 




number of par- 
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X has units of length. Also, since the thickness of the sample (height 
in z-direction) is much smaller than its size L, B curr differs very much 
from the volume magnetisation (demagnetisation effect). So one have 
to estimate the field -B CU rr and compare it to B cxt . The most strin- 
gent limitation is set by the requirement, that the largest possible 
value the local field -B CU rr can assume should be much smaller than 
.Bcxt. This yields 9 \x\k F ln(k F L) < 1 , or \x/xl\ < \n(k F L)/k F r e , 
where r e = e 2 /mc 2 is the classical electron radius. (Though, in some 
cases, e.g. in a square geometry, the condition is less severe, namely 
\x/Xl\ <S L\n(k F L)/r e .) As an example consider a disc geometry, 
where at paramagnetic peaks x/\Xl\ — (k F L) 2 A/T. If \ F ~ 1CT 6 cm 
and k F L ~ 100, then, as long asT> 10~ 3 A, the self-consistent treat- 
ment is not necessary. If the condition is not satisfied, one cannot 
assume anymore a given external field, but should solve the entire prob- 
lem self-consistently. A self-consistent treatment leads to such possibil- 
ities as Meissner effect and orbital ferromagnetism, both in bulk sam- 
ples and rings . 2 ' 3,19 ~ 22 An interesting possibility appears to exist in the 
strip geometry, considered in Section 3. Here the linear susceptibility, 
Eq. (18), is Xo — \XL\\Jk F L y . However, the differential magnetic sus- 
ceptibility, Xd, in a finite field B < mcv F /eL y , turns out to be of order 
\XL\(k F L y ) 3 / 2 (see Ref. 12, where a strip with harmonic confinement 
was analysed). 23 If |xo| <C \ F /4:ir\n(k F L) but \xd\ > \ F /Anln(k F L), 
then one encounters a situation similar to the one which can occur in 
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the de Haas-van- Alphen effect and leads to a break up of the sample 
into magnetic domains. 24 
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Figures 

Fig. 1. Linear magnetic susceptibility in a disc geometry, calculated at 
T = 0.1 A and normalised to \xL\(k F R) 2 . 

Fig. 2. Linear magnetic susceptibility in a disc geometry, calculated at 
T = 5 A and normalised to |xl|(^f-R) 3//2 - Analytic result obtained from 
Eq. (26) (solid line) is compared to numerical one based on Eq. (35) (dots). 

Fig. 3. Linear magnetic susceptibility in a square geometry, calculated at 
T = 0.1 A and normalised to \xL\(k F R) 3/2 . 

Fig. 4. Linear magnetic susceptibility in a square geometry, calculated 
at T = 5 A and normalised to \xl\- Analytic result obtained from Eq. (28) 
(solid line) is compared to numerical one based on Eq. (35) (dots). 
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1 (by E. Gurevich and B. Shapiro) 




Fig. 2 (by E. Gurevich and B. Shapiro) 



